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'^ ' Abstract 

- It has recently been shown that through the Tiy decay mode the heavy charged 

rS I Higgs bosons H can be discovered at LHC over a vast region of the parameter 

C^ '_ space previously thought unaccessible. If the charged Higgs bosons are discovered 

in this decay mode and the measurements of its mass and tan f3 are made, it allows 
us to make greater use of the dominant but QCD background contaminated bt decay 
mode by cuts suppressing phase space for the background. In this paper we study 
the top quark polarisation in the bt decay mode. When tan/? is either moderately 
large (tan /3 > 9) or small (tan /3 < 4) , the correlation effects can be used as a signal 
of H^ production. We outline the analysis procedure for various decay signatures. 
As an example, we demonstrate the application of our procedure to the dominant 
production process gb —>■ tH~ at LHC. The effect is clearly observable for H^ mass 
between 200 GeV and 500 GeV in both leptonic and hadronic channels, but at 500 
GeV the boost on the top quark starts to prevent the reconstruction of invariant 
masses, and hence the observation of the correlation effect, in the hadronic decay 
channel. 



1 Introduction 

The charged Higgs bosons H^ are a central ingredient of the Two Higgs Doublet Model 
(2HDM) which is not only an attractive alternative to the Standard Model single Higgs 
doublet but is also inherent in the supersymmetric models, notably Minimal Supersym- 
metric Standard Model (Mssm). Within the framework of MsSM, recent experimental 
analyses [|I], point to the heavy {Mh± ^ rrit) mass region, where they can not be pro- 
duced in top quark decay, and their dominant decay if^ -^ ht suffers from large QCD 
background. 

In the author has shown that the tv decay mode of the heavy charged Higgs 
bosons H"^ offers a striking signature of high p^ tau and large missing energy through 
which if ^ can be discovered readily at LHC, provided that the top quark momentum in 
the dominant production process gh — >■ tH~ is constrained. Although the signal cross 
section falls markedly at large Mh± and small tan/?, it was claimed that for tan/3 > 3 
and Mh± < 1 TeV the signal is sizeable enough after selection cuts to compete with the 
irreducible background and the finite LHC luminosity. 

Given now that heavy H^ can be discovered in the tu decay mode for much of the 
parameter space, and that Mfj± and tan/3 can be measured to certain extent using the tau 
and missing transverse momenta distributions and the total cross section, it is desirable 
to be able to narrow down the parameters further by utilising the dominant decay mode 
ht. 

The ht decay mode of if^, following production from the parton level process gh -^ 
tH~ , has been studied extensively in the literature |Q from the viewpoint of the search 
for heavy H^. The utility of the ht decay mode is limited for this purpose because of 
the contamination from QcD background. However, given that Mh± and tan/3 can be 
measured in the tu decay mode, this can be used to constrain the final state. 

In this paper we study the effect of top quark polarisation from H"^ decay on parton 
level observables, noting that the top quark is polarised by the chiral couphng of if^ 
to ht, for either large or small tan/?. To be more precise, the two chiral components of 
the H^ coupling to ht are proportional to rrif, tan (3 and rrit cot /3 and the two balance 
out when tan/3 = ^Jm^Jmy, ~ 6.4 {m\^ = 4.25 GeV). Far away from this central value, 
tan^/3 <C rri^/ml or tan"^/? 3> m^/ml, the couplings are dominantly left or right leading 
to characteristic correlation patterns. 

Our discussions are applicable to any class of production processes, but as an example 
for studying the effect of boosts on H"^ momenta we adopt the gh —>■ tH~ production 
process at LHC. 

When tan/3 ~ rrit/mb, the production cross section in the gh -^ tH~ reaches the min- 
imum, while the correlation pattern will resemble the background, and our analysis will 
not be of much use. If the background reduction is almost perfect before the polarisation 
analysis, the polarisation measurement can be used to measure tan/3 accurately. This is 
presumably unlikely. 

The use of right-handed tau polarisation as a characteristic signal of H^ decay was 
studied in |^. In addition, we note that if the tau is produced in the decay of the top 
quark through H"^, the distribution of invariant mass squared Mn^^ is flat since there 
is no momentum correlation between the h from the decay of the top quark and the 
r. This does not apply to the background from W^ decay, and if we know the exact 
distribution, we can use it to measure the signal to background ratio and thus confirm 



the tau polarisation analysis. 

Our current work is in a similar vein. As in ||^, the technique is not expected to be 
of great use in reducing the background, but once an excess of signal over background is 
observed near if ^ mass, the momentum correlations resulting from top quark polarisa- 
tions can be used as a check for the signal presence in the sample and its proportion to 
the background. 

We outline the analysis procedure for various decay signatures. 

Finally, we note that a preliminary study of the top quark polarisation in H^ decay 
was presented in [^ . Our present analysis is different in several ways. Firstly, the primary 
purpose of is in measuring tan/3, rather than in distinguishing the signal from the 
background. Thus their analysis is applicable in a range of tan (3 complementary to the 
range emphasised here. Secondly, the variable utilised in 0, namely the leptonic energy 
in the H"^ rest frame, is not well defined unless H^ is produced at an e^e~ collider 
and decays leptonically with no other source of missing energy. In any case, for the 
intermediate region of tan /3 the production rate at hadronic colliders is small. 

2 Calculation 

We adopt Mh± and tan (3 as the two Higgs sector parameters. When calculating the width 
of H^ in order to simulate the finite width effect, we assume that the ht branching ratio 
is 100%. The top branching ratio to hW^ is also taken to be 100%, and W^ branching 
ratio to a pair of fermions //' is divided into quarks and leptons by the tree level ratios 
of 1/3 for each quark pair and 1/9 for each lepton pair. Note that unlike in reference 0, 
'lepton', £, always includes tau. 

We restrict our discussions to the tree level. The dominant part of the radiative 
corrections [0, which are known to be significant, can be absorbed in terms of the running 
mass of the bottom quark and therefore does not affect the shapes of our distributions. 
We simulate the jet showering by imposing cuts between parton directions for the case 
of the hadronic decay of the top quark. 

For simplicity we take the top quark mass to be 175 GeV for both Yukawa coupling 
and the kinematics, and the bottom quark mass to be 4.25 GeV, again for both Yukawa 
coupling and the kinematics. The fermions from W^ decay are considered massless. 

The electroweak parameters are aEM = 1/128, siv? 6w = 0.2315, Mz = 91.187 
GeV, Mw = Mz cos 6w ~ 79.94 GeV. The Cabbibo-Kobayashi-Maskawa matrix ele- 
ment V(?j^jy[[6t] is taken to be 1. 

The tree level matrix elements squared are simple and are calculated by hand. Con- 
sistency checks are carried out between different limits. The numerical integrations are 
carried out using VEGAS 0] and the total decay width, in the narrow width limit of the 
propagators, agrees with the width calculated by hand. 

3 Analysis 

The decay matrix element squared, corresponding to the Feynman graph of figure |l|, is 
given by: 
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[diM^i -ml- pf) - 2gbgtmbmt] + —^ — -{gt^t - dPt) 

■^Pt ■ Pi 

We have defined gt = {rrit cot (3/ Mw) and g^ = {mbta.nP/M]v). Qt and Qw^ are Breit- 
Wigner propagators for the top quark and W^ , respectively. In our analysis we also 
generate the H^ mass with a Breit-Wigner distribution. This formula essentially carries 
all the information about momentum correlations. The effect of top quark polarisation is 
manifested in the second term in large square brackets which is proportional to {2pi -pf). 
If the two chiral couplings were equal, the term vanishes in the limit of on-shell top quark. 
In the limit of on-shell top quark and W^ boson, the matrix element squared, and 
hence the differential decay rate, depend only on two independent dot products (or in- 
variant masses), which we can take to be {2pi -pi) and {2pb-p() since these two quantities 
are measurable, as invariant masses, in the leptonic decay of the top quark, if bottom 
quark charge tagging is available. We have: 

2pb ■ p^ = m] -ml- M^± - 2pb ■ pi 

2pfpe = M^i+2pb-pi (2) 

and hence the above statement holds. Furthermore, from equations (HD and (^ we see 
that in the limit of non-chiral coupling, the dependence on {2pi ■ pi) vanishes altogether, 
as would a general non-chiral background. Hence the dependence of the differential decay 
rate on {2pi ■ pi) would be an indication of the presence of if ^. This would not hold if 
the two coupling strengths were similar, g^ ~ g^. This implies tan^/3 ~ mf/ml. 

In experimental analyses, a possible approach would be to fit the multi-dimensional 
differential decay distributions with theoretical expectations for signal and background. 
However, for our purposes, let us define the following dimensionless quantity in the limit 
of on-shell W^: 

\= \ ^[5^] ~ "^fe ^ ), 2pfc ■ pe ^ 2pi ■ pi 

'Ml,^-ml + Ml, ^'2pb-p, + Ml, °2p,.p/ ^'^ 

From (|I]), the differential distribution of A translates directly to the chiral coupling 
strengths. Aq is the normalisation which is a function of the masses. Given that Mu± 
needs to be known approximately for our analysis to be of much use, we can define Aq 
for instance as follows: 

^ ^ mj - M^± _ (2pb-p£)max ,^s 

M|^± - ml ^ {2pi ■ p^)max ' 

The dependence on A observed will be the sum of the signal and the background which, 
after the signal selection cuts, is compatible with H^ production. If the background does 
not involve chiral interactions, the distribution is fiat, and even when the background does 
involve chiral interactions, the distribution can be analysed as the linear combination of 
the signal and background expectations. One can compare the ratio thus obtained with 
the signal to background ratio observed in order to confirm the presence of H^ signal. 

In the next chapter we study the exact definitions of A for different final states. 

4 Results 



4.1 Leptonic decay with charge tagging 

The simplest case is the leptonic decay of the top quark, with the charge tagging of both 
bottom jets and the lepton. This allows us to assign all momenta correctly. The size of 
the sample will be limited, since in addition to the leptonic branching ratio of the top 
quark, the rate is reduced further by the requirement of bottom quark charge tagging, 
the efficiency of which one might roughly estimate to be around (10%)^ = 1%. 

In figure |^ we show the differential distribution of A at M^± = 500 GeV. As stated 
in the previous chapter, the distribution is flat for tan/5 = ^mt/rrii, ^ 6.4, and rapidly 
approaches the large and small tan/3 limits away from this central value. Let us concen- 
trate on the two representative extreme values tan/? = 1.5 and 30, and the central value 
tan/3 = 6.417 from here on. 

The flnite energy resolution of detectors can lead to a flnite spread in the invariant 
mass distribution [y]. The binning in our distributions of A A = 0.1 are fairly broad, 
corresponding to roughly 0.1 times the energy scale considered. For the quantity {2pt -pi) 
in the denominator of (^, this corresponds to the energy spread of O.lmt ~ 15 GeV. This 
is similar to the resolution obtainable at LHC ||T0[| . 



The dependence on Mh± is studied at tan/3 = 1.5 and tan/3 = 30 in flgure ^. We see 
that the scale normalisation of (|) is satisfactory, and there is little residual dependence 
on Mh±- 

4.2 Leptonic decay without charge tagging 

Once we drop the bottom quark charge tagging requirement, we need to consider a new 
definition of b and b in the procedure outlined above. It is possible that the signal selection 
procedure has already picked out one or the other of the two bottom jets as belonging to 
the decay of the top quark. If this is the case, the assignment is again trivial. From here 
on we assume that this is not the case. 

If Mh± ^ rrit the procedure is straightforward, as b and the top quark decay products 
point back-to-back, even when we consider the boost from the production process. In 
this case 2ph ■ Pe ^ m^ — M^ and 2pi ■ pi = 0{Mjj±) is much greater than the other dot 
product. The author has verified that the distributions look identical to figure 0a if we 
simply take {2pb ■ pe) to be the smaller of the two dot products. 

For general Mh±, particularly when (M|f± — m^) ~ {mf — M^±) at Mh± ~ 250 GeV, 
the two dot products have similar distributions. In this case the charge tagging discussed 
above is an effective strategy. Here we proceed by summing over the two distributions 
corresponding to the 'right' and 'wrong' assignments of bottom quark jets. 

The resulting distributions are shown in figure ^ for Mh± = 200 GeV. This is already 
sufficiently below the contaminated region. Note that the normalisation is to two. Com- 
paring with figure 0a we see that the region A < 1 is similar in two distributions, but 
figure ^ has tails at A > 1 where most of the wrong assignment lies. 

This result indicates that the mis- assignment of bottom quarks, for example from 
production processes, is unlikely to have significant effect on our analysis, because of the 
different kinematic limits on the invariant masses. 

We see that the flatness of the 'background' tan/? = ^yrrit/mb case is retained here. 



4.3 Hadronic decay 

If the top quark decays hadronically, the mass reconstruction allows one to assign b and 
b correctly to the two bottom quark jets. We assume that this has been done already in 
the signal selection process. We also assume that bottom quark jets are distinguishable 
from light quark jets. However, the i and z/ of equation (^ are now indistinguishable 
unless we tag the charges of light quark jets as well as the bottom quark jets. 
Here we choose to proceed by defining i to be the jet which minimises pb ■ Pjet- 
The result is shown in figure ^ The distribution is somewhat flatter than the distri- 
butions in the leptonic channel, figures H and §, reflecting the difficulty with which the 
light quark jets can be assigned to the parent quarks. Mh± = 500 GeV is adopted, but 
the distribution for different Mh± are found to be similar. 

The problem of assigning the jets correctly can be seen in figure |^ which compares 
the differential distributions in invariant masses M?^^-, and Mn^-,y The distributions are 
roughly independent of either Mh± or tan (3. We see that the two dot products have simi- 
lar profiles, and our procedure can not escape contamination from the wrong assignment. 

4.4 The effect of cuts for the hadronic decay 

The above discussions for the hadronic decay assume that the parton directions are known 
from the jet distributions. In practice, the jets can overlap when either top quark decay 
products being produced from relatively light H^ interfere with the b jet, or when the 
top quark is boosted in heavy H"^ decay and the decay products are clustered near each 
other. These cases lead to ambiguous definitions of jet momenta. 

Let us study the effect of jet showering by introducing cuts on jet separations. 

A typical azimuthal-pseudorapidity separation A_R = A/(A0p~-i-~(^^?)^ required to 
define jets and reduce the inter-jet interference might be taken to be 0.7. In terms of the 
azimuthal angle at zero rapidity this corresponds to 0.7 radians, or 40 degrees. 

Let us proceed by replacing AR by the angle between two jets in the rest frame of 

We do not consider the effect of nonzero H^ momenta. This is sensitive to the selection 
cuts imposed in order to reduce the background. 

The ratio of events passing the cuts is shown in figure ^ The suppression is severe 
at large Mh± where the production cross section also falls. From the events which do 
pass the separation cuts, it is difficult to extract correlation effects. At Mh± = 1 TeV, 
for example, the effects vanish altogether. 

The distribution after cuts is shown in figures ^ for Mh± = 200 GeV and 500 GeV, 
to be compared with figure ^ before cuts. The normalisations are to unity before cuts. 

The effect of cuts is found to be significant, but not prohibitive in studying top quark 
polarisation. 

4.5 Production processes 

In order to study the effect of boosts on H^ momenta from the production process, we 
consider the production process gb -^ tH^ . The distributions are expected to be modified 
slightly because of the boosts on H^ momenta acquired from the production. For the 
calculation of the production process we proceed as in B. 



In figures || we show distributions analogous to ^ for the hadronic decay after cuts. 
The cuts are now in terms of AR in the colhder frame rather than the angle between 
jets. The normalisation is to unity after cuts for each distribution. We see that the 
distributions are not significantly affected by the boost. 

The detector cuts on pseudorapidity, and any other separation cuts involving jets 
arising from the production process, for example the top quark in gb —>■ tH^, can only 
affect the normalisation of the distribution. The separation cuts between the H^ decay 
products are significant only because they directly affect the range of invariant masses 
accessible for our analysis. 

For the sake of completeness, the author has run an analogous simulation for the 
leptonic top quark decay with the cuts AR > 0.7 between the bottom quark jets and 
AR > 0.4 between the lepton and the bottom quark jets. The effect has been found to 
be small, as expected. 

5 Conclusion 

In this paper we studied the effect on top quark polarisation of the chiral charged Higgs 
boson H^ coupling to fermions. We discussed the optimum variables A which can be 
adopted in both leptonic and hadronic decays of the top quark. 

The variable A we adopted assumes that Mh± is already approximately known by the 
time we apply this analysis, either by the observation of a peak in invariant mass of tb, 
or by the analysis of H^ decay to tu as mentioned in [^. 

For the leptonic decay of the top quark, A can be defined such that the distributions 
have the desired features of a clear momentum correlation effect and fiat distributions 
for the background. Charge tagging the bottom quark jets and the lepton helps when 
Mh± ~ 250 GeV. 

The correlation effect is significant both for large and small Mh±. 

For the hadronic decay of the top quark, the momentum correlation effects are less 
pronounced, but it is still strong enough to be used as a signal for H^ production. 

The distributions are affected significantly in the hadronic decay case by the require- 
ment of jet separations. This is especially the case for heavy H^. 

The boost on if^ from the production process is found to give only a small effect, 
in the case of gb —>■ tH^ production mode. We expect our analysis to be general to all 
production modes, at LHC and beyond. 

In conclusion, top quark polarisation is useful for confirming the presence of H^ in a 
sample of final states containing tb in experimental analyses of near future, at tan /3 ^ 9 
and tan/5 < 4, and up to Mh± ~ 500 in both hadronic and leptonic channels. 
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Figure captions 



Figure 1: The tree level Feynman graph for the decay H^ -^ tb and subsequent top quark 
and W^ boson decays. The H~ decay can be obtained trivially from H^ decay. 

Figure 2: The differential distribution of A for the leptonic decay of the top quark, with 
charge tagging. Normalisation is to unity. Lines were used instead of histograms for the 
sake of clarity. The binning width is 0.1. (a - top) Mh± = 500 GeV, for five different 
values of tan/3, (b - bottom) tan/? = 1.5 and tan/? = 30, for three different values of 
Mh±. The tan/? = 30 distributions have negative gradients. 

Figure 3: The differential distribution of A for the leptonic decay of the top quark, without 
charge tagging. Mh± = 200 GeV, for three different values of tan /3. Normalisation is to 
two. 

Figure 4: The differential distribution of A for the hadronic decay of the top quark, before 
cuts. Mfj± = 500 GeV, for three different values of tan /?. 

Figure 5: The differential distribution of two invariant masses M?,i and M?i^-n, corre- 
sponding to 'b-nu' and 'b-lep' respectively in the figure legends, in the hadronic decay 
of the top quark. M//± = 500 GeV, tan/? = 1.5. 

Figure 6: The ratio of H"^ decay events which pass the jet separation requirement, for 
three different values of tan/3 and 200 GeV < M^± < 1 TeV. 

Figure 7: The differential distribution of A for the hadronic decay of the top quark, after 
cuts, (a - top) Mjj± = 200 GeV; (b - bottom) Mh± = 500 GeV. The normahsation is to 
unity before cuts. 

Figure 8: The differential distribution of A for the hadronic decay of the top quark, after 
cuts, convoluted with gb -^ iH^ production process, (a - top) Mh± = 200 GeV; (b - 
bottom) Mh± = 500 GeV. The normalisation is to unity after cuts. 
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